Selection of special orientations in relativistic collisions of deformed heavy nuclei 
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We have studied U+U collisions at ^/sjvjv = 200 GeV using Monte Carlo Glauber, UrQMD and 
AMPT models. We find that it is possible to separate central tip-tip events as well as central 
body-body events on the basis of cuts on multiplicity and magnitude of the reduced flow vector. 



PACS numbers: 25.75.-q, 24.10.Lx 



I. INTRODUCTION 



One of the recent surprises about the properties of the 
QCD matter formed in Au+Au collisions at RHIC is that 
it appears to flow with near-zero viscosity: an almost per- 
fect liquid [EH&HHH. ThiS 

is in contrast with prior 
expectations of non-interacting gas-like behavior, and is 
based in part on approximate agreement between results 
from non- viscous hydrodynamic calculations and experi- 
mental data on elliptic flow. Azimuthal anisotropy in the 
transverse plane, or elliptic flow (^2), is sensitive to the 
early pressure developed in the system and its equation 
of state. There is a close connection between elliptic flow 
and the initial spatial eccentricity (e) of the transverse 
overlap region of the two colliding nuclei. Non- viscous 
hydrodynamic calculations predict constant vijt over a 
wide range of impact parameters While the data 
remain well below the non-viscous hydrodynamic values 
at lower beam energies, central collisions of Au+Au at 
\J s nn = 200 GeV just reach the prediction 0,0]. 

In relativistic uranium-uranium (U+U) collisions, 
there is the potential to produce more extreme condi- 
tions of excited matter (higher density and/or a greater 
volume of highly excited matter) than is possible using 
spherical nuclei like Au+Au or Pb+Pb at the same in- 
cident energy [1, H, [HI, E3 • However, this potential 
can be partly lost if we have limited capability to dis- 
tinguish experimentally between different collision orien- 
tations when the ions interact near zero impact param- 
eter. The prior studies cited above have devoted atten- 
tion to this problem, and while there is agreement that 
U+U collisions offer worthwhile advantages over Au+Au, 
quantitative particulars have yet to be worked out in de- 
tail [n], HH, EH- Furthermore, U+U collisions offer the 
opportunity to explore a different and wider range of ini- 
tial eccentricities than is possible in the case of Au+Au. 
Among the possible collision orientations, of particular 
interest are the "tip-tip" orientation, in which the long 
axes of both deformed nuclei are aligned with the beam 
axis, and the "body-body" orientation, in which the long 
axes are both perpendicular to the beam axis and par- 
allel to each other. In this paper, we report a study 
of the selection of these orientations based on a Monte- 
Carlo Glauber model [HI, [HI, Ultra-relativistic Quan- 
tum Molecular Dynamics (UrQMD 1.3) flUl. and A Multi 



Phase Transport Model (AMPT vl.ll-v2.11) [L($. 

The Glauber model is solely based on collision geome- 
try [li|. The UrQMD model is a microscopic transport 
theory based on the covariant propagation of all hadrons 
on classical trajectories in combination with stochastic 
binary scatterings, color string formation and resonance 
decay [la ]. We use this model with default settings. 
AMPT is a hybrid model. It uses the heavy-ion jet inter- 
action generator (HIJING) to generate initial conditions, 
Zhang's parton cascade (ZPC) for the partonic scatter- 
ing and hadronization, and a relativistic transport (ART) 
model for hadronic interactions and freeze-out [161 ]. We 
use this model with string melting, where all excited 
strings that are not projectile and target nucleons with- 
out any interactions are converted to partons according 
to the flavor and spin structures of their valence quarks. 
This option with appropriate partonic cross-section gives 
higher value of elliptic flow, close to the experimental 
value, than without string melting [l(|. We keep all other 
options at their default settings. 

The article is organized as follows: in Sec. HI] we briefly 
review the relevant quantities and outline our calculation, 
and in Sec. IIII1 we present the selection procedure for 
tip-tip and body-body events and discuss the results. In 
Sec. IIV1 wc summarize our findings. 



II. CALCULATION AL FRAMEWORK 

We represent the quadrupole deformation of the 
ground-state uranium nucleus in the standard [ItJ way: 
we take a Saxon- Woods density distribution with sur- 
face thickness a = 0.535 fm and with R = i? sp (0.91 + 
0.27 cos 2 9), where 9 is the polar angle relative to the long 
axis of the nucleus and R sp = 1.12 A 1 / 3 - 0.86 A" 1 / 3 
fm [l?]]- The small hexadecapole moment of the ura- 
nium nucleus is neglected. This yields i?iong/-Rpcrp ~ 
1.29 [l2j|- The orientation of the first and second nucleus 
in the colliding pair is fixed by the two angles (9 p ,cj> p ) 
and (9 t ,4>t), respectively. The angles 9 p and 9 t describe 
the orientation of the long axis relative to the beam di- 
rection, and they are uniformly distributed in [0, 7r/2]. 
The azimuthal angles <f> p and 4>t describe rotations about 
the beam direction, and they are uniformly distributed 
in [0, 27r]. The Monte-Carlo Glauber simulation is as 
described in Rcf. and chareed multiolicitv in the 
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central pseudo-rapidity region is parametrized using the The area of the transverse overlap region is given by 
approach of Ref . 18| : 



dN, 



ch 



drj 
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(1) 



where n pp = 2.19 and x = 0.15 are obtained by fitting 
the PHOBOS data [U] for Au+Au at 200 GeV [3. 
Here, N^ and N w are the number of binary collisions, 
and number of participant (wounded) nucleons, respec- 
tively. The multiplicity is then distributed according 
to a negative binomial 20] in order to get a realistic 
distribution. To convert between track densities in ra- 
pidity y, and pseudo-rapidity r/, we use the approxima- 
tion dN ch /dy w 1.15 dN ch /drj [2l[ in the Monte-Carlo 
Glauber calculations. The UrQMD and AMPT codes 
have been modified to simulate the deformed uranium 
nuclei. 

In this paper, "ideal tip-tip" and "ideal body-body" 
refer to the configurations in which the long axes of both 
deformed nuclei are aligned with the beam axis at zero 
impact parameter (9 p = 6t = 4>p = 4>t = 0, 6 = fm) 
and in which the long axes are both perpendicular to 
the beam axis and parallel to each other at zero impact 
parameter (fL = Q t = n/2, (f> p = <f) t = 0, 6 = fm), 
respectively [l, [l2[ • As the statistical weight of the above 
mentioned ideal configurations is negligible, "tip-tip" and 
"body-body" are used to refer to the configurations close 
to "ideal tip-tip" and "ideal body-body" , respectively, 
for selection purposes. 

In the standard way of calculating the spatial eccen- 
tricity of the transverse overlap region, e stc i, the minor 
axis of the ellipse representing the transverse overlap re- 
gion is taken along the direction of the impact parameter. 
However, the minor axis may not be along the direction 
of the impact parameter due to fluctuation in the par- 
ticipating nucleon positions [22| . The spatial eccentricity 
calculated by taking into account the rotation of the mi- 
nor axis is referred to as the participant eccentricity, e par t, 
and given by Ref. [22J: 



tpart 



4er 2 

xy 



(2) 



with 



ol = (x 2 ) (x) 2 



°i = (y 2 ) (y) 2 , (3) 

(?xy = (xy) - {x} (y) , (4) 

where x and y are the transverse coordinates of the par- 
ticipant nucleons and angle brackets denote averaging 
over participant nucleons in a single event. For spher- 
ical nuclei, the difference between e stc i and e par t is sig- 
nificant only in smaller systems or smaller overlap re- 
gions [22| . However, in case of collisions of deformed 
nuclei, the transverse overlap region may not be an el- 
lipse depending on the collision configurations. Because 
of this, the difference between e st d and e par t is significant 
even in central collisions. We use e Da rt in our calculations. 



S = tt \ I a 2 (Jy 



(5) 



The event flow vector for the second harmonic, Q 2 , is 
defined as 1231 



Q 2 = ^2 iw k cos 2cj) k + ^2 jw k sin 2<j) k . (6) 



We also use the corresponding reduced flow vector, whose 
magnitude is 
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where <\> k is the azimuthal angle of a particle, w k is a 
weighting factor for that particle, and the summation is 
over all n cri charged particles in the event. The magni- 
tude of the reduced flow vector, q 2 as defined above, is 
used in this work to establish a measure of flow strength 
where the trivial dependence on n cri is removed [2l| . The 
transverse- momentum- integrated elliptic flow. is pa- 
rameterized as v 2 = 0.034 e part (dN^/dy) 1 / 3 [2J], where 
dN^/dy is the central rapidity density of the charged 
particles. This expression is known to give a good de- 
scription of Au+Au collisions at ^snn = 200 GeV and 
we assume here that it also holds for U+U collisions. 

To include v 2 in the Glauber model, azimuthal angles 
cf> for 2(dN c h/drj) tracks are generated according to the 
above paramctrization, and we take Wk = 1. In the 
UrQMD and AMPT models, w k = for < 2.0 
GeV/c and w k = 2.0 GeV/c otherwise, where p-j; is the 
transverse momentum of the fc th particle in the event. 

The quantities S and e part for all the above- 
mentioned models were calculated on the basis of the 
Glauber model. The transverse particle densities, 
(I /S) (dNc h /dy) [HG1 for ideal cases are: 



1 dNch 
S dy 



43.9 fm ideal tip-tip , 

31.0 fm~ 2 ideal body-body . 



(8) 



III. SELECTION PROCEDURE 

The calculations based on the models show that the 
U+U central tip-tip configuration has the smallest over- 
lap area, similar to central Au+Au, and results in the 
highest multiplicity. It is a good candidate to create the 
highest possible (1/ ' S)(dN c h/dy) at any given beam en- 
ergy. The higher the multiplicity, the better the appli- 
cability of hydrodynamics calculations [25[ . The central 
body-body configuration also yields higher multiplicity 
than average, but smaller (1/ S)(dN c h/dy), similar to cen- 
tral Au+Au, due to its large overlap area. However, it 
produces larger elliptic flow at the same (l/S)(dN c h/dy). 
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The tip-tip configuration yields higher (l/S)(dN c \ i /dy) at 
similar e, and the body-body configuration yields higher 
e at similar (l/S)(dN c i 1 /dy), compared to Au+Au at the 
same beam energy. For comparison of the selected events 
with the ideal cases, we fitted the S and (1/ ' S)(dN c h/ dy) 
distributions for ideal tip-tip and ideal body-body config- 
urations with Gaussians. The quantity q 2 for ideal tip-tip 
is fitted according to [2l| : 



dP 
q 2 dq 2 

with 



f 



cxp 



v 2 2M 



4 



2*1 



<72«2 



M 



(9) 



^ = 0.5(1+52), (10) 

where Jo is the modified Bcssel function, M = dN c h/dr], 
and v 2 is the elliptic flow. The quantities v 2 and g 2 are 
taken as free parameters in fitting. The participant ec- 
centricity, e par t , is fitted similarly with different values of 
g 2 . For ideal body-body, a Gaussian fits well. 

The following characteristics of the tip-tip and body- 
body configurations are used for selection: 

{largest dN c h/drj (central collisions) , 
smallest q 2 . i.e., smallest e par t, and 
smallest S ; 

(11) 



and 



body-body 



largest dN c h/dr] (central collisions) 
largest q 2 , i.e., largest e part , and 
largest S . 



(12) 

The selected events are compared to the ideal cases. 
While 5* and e part cannot be directly measured, they are 
compared to the ideal cases for identification purposes. 



A. Tip-tip 

The quantity N w is approximately the same in tip- 
tip and body-body configurations because both have full 
overlap. However, TVb is larger in tip-tip than body-body. 
Because of the contribution from Ab, the multiplicity 
is larger in tip-tip than body-body collisions. Thus, it 
is possible to select tip-tip on the basis of multiplicity 
only. However, for a variety of reason, this approach has 
relatively poor efficiency, and instead we have used both 
multiplicity and the reduced flow vector. 

To select tip-tip events, first we take the top 1% 
dN c h/dri to select full overlap configurations. These are 
already rich in tip-tip, but a further selection of the bot- 
tom 25% q 2 is taken from the above subset to further en- 
hance the tip-tip configuration. Figs.[TJ[2l and[3]show the 
event distributions for the above three models, after cuts, 
compared with the ideal tip-tip case. The transverse par- 
ticle density, (1/ S)(dN c \ l /dy), obtained after these cuts 
is 40 fm _ , which is significantly larger than Au+Au 
at the same beam energy. 



Event distribution for tip-tip selection: Glauber 
(a) 




FIG. 1: (Color online) Distribution of (a) transverse over- 
lap area, (b) eccentricity of the transverse overlap region, (c) 
magnitude of the reduced flow vector, and (d) transverse par- 
ticle density from Monte-Carlo Glauber. The y-axes have 
arbitrary units. The dotted and dashed curves are fitted for 
ideal tip-tip and ideal body-body configurations respectively, 
compared here with the selected tip-tip events (solid lines). 



Event distribution for tip-tip selection: UrQMD 
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FIG. 2: (Color online) Same as Fig.[l]for UrQMD. 



B. Body-body 

Here, we first take the top 3% dN^/dr) to select full 
overlap configurations. Then we take the top 1% q 2 to 
select body-body events. The (1/ S)(dN c \ l / dy) obtained 
after all these cuts is ~ 33 fm -2 . Figs. |4j [5j and [6] show 
the event distributions, after cuts, compared with the 
ideal body-body case, for the three models. 

We have also tested an alternative procedure: take the 
top 3% dNch/dri to select full overlap events. Then, take 
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Event distribution for tip-tip selection: AMPT 
(a) 
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FIG. 3: (Color online) Same as Fig. Q] for AMPT. 
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Event distribution for body-body selection: UrQMD 
(a) 2sF 
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FIG. 5: (Color online) Same as Fig.Q]for selected body-body 
events in the UrQMD model. 



the lower 50% dN c h/dr] from the above subset to select 
mostly body-body events. This sample is further purified 
by taking the top 5% q 2 . However, this method is less 
efficient than the one discussed above. 



Event distribution for body-body selection: Glauber 
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(Color online) Same as Fig.[T]for selected body-body 
events in the Glauber model. 



C. Discussion 

The center of the S distribution in the ideal tip- 
tip and the ideal body-body from AMPT (Fig. EJJa) 
and Fig. Ela)) differs from that for the Glauber model 
(Fig. (Ha) and Fig. HJa)) and for UrQMD (Fig. H» 
and Fig. 0a)) because of a different parametrization of 
the nuclear radius used in AMPT compared with the 



Event distribution for body-body selection: AMPT 
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FIG. 6: (Color online) Same as Fig.[T]for selected body-body 
events in the AMPT model. 



other two models. The e part distributions for the selected 
events from all models are in good agreement with the 
ideal cases for both tip-tip (Fig.EEfb), Fig. 12(b), Fig.^b)) 
and body-body (Fig.Hfb), Fig.[5(b), Fig.|5(b)) selections. 
The quantity q2 is sensitive to e par t and because of this, 
the hard cuts on q 2 bring the selected sample closer to 
the ideal values. 

The dN c h/dr] cut selects mostly tip-tip events, with 
some body-body as well as other event categories having 
smaller e par t but larger S, depending on the orientation 
of the nuclei. The additional cut in q 2 helps to remove 
most of the body-body events, but still leaves some events 
with smaller e par t but larger 5". That is why the S distri- 
bution for the selected events is a little different from the 
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ideal cases. The selection purity can be further improved 
by applying an even harder cut in dN c h/dr]. A similar 
argument is also applicable to the body-body selection. 

IV. SUMMARY AND CONCLUSION 

U+U collisions at maximum RHIC energy are pre- 
dicted by three rather different models to produce more 
extreme excitation of nuclear matter than has been 
achieved to date using Au+Au. The expected benefits 
of colliding U+U include a closer approach to the condi- 
tions where hydrodynamics should be applicable, and the 
opportunity to explore a different and wider range of ini- 
tial eccentricities than is possible in the case of Au+Au. 
However, we can achieve the full potential benefits of 
U+U collisions only if we are able to select favored colli- 
sion orientations. 

The parameter space explored by colliding deformed 
nuclei is more varied than for spherical nuclei. The tip- 
tip U+U configuration offers higher (l/S)(dN c h/dy), and 
lower V2 (similar to central Au+Au). The body-body 



configuration gives higher t>2, and lower (l/S)(dN c \ l /dy) 
(similar to central Au+Au). 

We have presented a prescription for selection of tip-tip 
and body-body configurations, based on the experimen- 
tal observables dN c h/drj and <j2- The e pa rt values for the 
selected events are in very good agreement with the ideal 
cases in all models, although S is somewhat different. In 
our previous study , we have estimated the benefits of 
U+U collisions including the effect of detector resolution. 
Here we have used a wider set of kinematic variables to 
demonstrate that the desired configurations can be se- 
lected efficiently. With the development of the Electron 
Beam Ion Source (EBIS) [26[, we expect U+U collisions 
at RHIC in the not-too-distant future. 
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